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The branching fractions of the decays B° — > K*°j and B + —+ K* + j are measured using a sample 
of 88 x 10 BB events collected with the BABAR detector at the PEP-II asymmetric-energy e + e~ 
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collider. We find B(B° -> K*°j) = (3.92 ± 0.20(stat.) ± 0.24(syst.)) x 10" 5 , B(B + -> ^'+7) = 
(3.87±0.28(stat.) ±0.26(syst.)) x 10~ . Our measurements also constrain the direct CP asymmetry 
to be -0.074 < A(B -> A* 7) < 0.049 and the isospin asymmetry to be -0.046 < A - < 0.146, 
both at the 90% confidence level. 

PACS numbers: 13.25.Hw,12.15.Hh,11.30.Er 



Within the Standard Model (SM), the decays 
B — ► A'* 7 proceed dominantly through one- loop b — > sj 
electromagnetic "penguin" transitions lj. Non-SM vir- 
tual particles may be present in these loops, changing the 
decay rates from the SM predictions. Theoretical calcula- 
tions of exclusive B — ► A*7 decay rates have large uncer- 
tainties due to nonperturbative hadronic effects 0, EJ El , 
limiting their usefulness for probing new physics. Pre- 
vious measurements 0, 0, °f the branching fractions 
are already more precise than SM-based theoretical esti- 
mates, and are in reasonable agreement with them. Cal- 
culations @> of the form factor for B — ► K*j can be 
tested using improved measurements of these branching 
fractions. 

Much of the theoretical uncertainty in the branching 
fractions cancels in the ratios defining the isospin asym- 
metry Aq- and the CP asymmetry A: 



t(b° -> r° 7 ) 

^0- - — =0 =To 



T(B~ -> K*~-f) 



T(B -> A 7) + T(B- K*~j) 



A = 



r(B -> a* 7 ) - r(B -> a* 7 ) 

T(B -> A* 7 ) + T(B -> A* 7 ) '' 



(1) 



(2) 



making them stringent tests of the SM. A further ad- 
vantage of these asymmetries is that some experimental 
systematic uncertainties cancel in the ratios. The SM 
predicts a positive value of A - between 5 and 10% [i"oj| . 
and \ A\ less than 1% ^lj. New ph ysic s contributions can 
modify these values significantly |lfl 1 1 11] . 

In this Letter, we present measurements of the exclu- 
sive branching fractions B(B° — » K*°j) and B(B + — > 
K* +r y), the isospin asymmetry (Ao_), and the CP asym- 
metries A(B° -> A*° 7 ) and -4(5+ -> iT* + -y). A"* refers 
to the A* (892) resonance throughout this paper. Inclu- 
sion of charge-conjugate decays is implied except in the 
definitions of A. This analysis uses (88 ± 1) x 10 6 BB 
events, from T(4S) decays, recorded by the BABAR de- 
tector 12]. An additional 10 ftr 1 of data, taken 40 MeV 
below the T(45) resonance, is used for studying non-i3 
continuum background. After B — ► A"*7 event recon- 
struction and background rejection, multi-dimensional 
extended maximum likelihood fits are used to extract the 
final results. 

We reconstruct B° -> A*° 7 in the A* -> A+tt", 
A°tt modes and B+ -> A*+ 7 in the A*+ -> A+tt , 
Ag7r + modes as described in detail in Ref. 00]. Re- 
constructed tracks are identified as final state tt^ and 



mesons by measuring the angle of the Cherenkov 
cone and energy loss along the track (dE/dx). The A^ 
candidates are composed from pairs of oppositely charged 
tracks with an invariant mass that is within 3.3er of the 
nominal A^ mass and with a vertex that is at least 0.3 cm 
away from the primary event vertex. The 7r°-candidate 
momentum vector is determined by a mass-constrained 
fit to pairs of photons, reconstructed from energy de- 
posits in the calorimeter that are not matched to tracks. 
The A and n candidates are combined to form A* can- 
didates, which are required to have invariant mass in the 
range 800 < M K7T < 1000MeV/c 2 . The primary-photon 
candidates are required to have high center-of-mass (CM) 
energy, between 1.5 and 3.5 GeV, and to satisfy addi- 
tional requirements designed to suppress the large ir° and 
•q background as described in Ref. |6( . 

The A-meson candidates are reconstructed by com- 
bining the A* and high-energy photon candidates. We 
define in the CM frame (denoted by asterisks) AA* = 
E* B — A^ cam , where E£ eam is the beam energy, known 
to high precision, and Eg = E* + E* K » is the energy 
of the A-meson candidate. We also define the beam- 

p's , where 



E* 2 

beam 



energy-substituted mass toes = 
p'* B is the momentum of the B candidate modified by 
scaling the photon energy to make A* + E* K , —El eam = 0. 
This procedure reduces the tail in the signal toes distri- 
bution, which results from the asymmetric calorimeter 
response. For signal decays, this "rescaled" tties peaks 
near 5.279 GeV/c 2 with a resolution of «3 MeV/c 2 and 
AA* peaks near MeV with a resolution of «50 MeV. 
We consider only candidates with toes > 5.20 GeV/ c 2 
and I AS* I < 0.3 GeV. 

Background events arise predominantly from random 
combinations of particles in qq production (q=u,d,s,c), 
with the high-energy photon originating from initial-state 
radiation or from tt° and 77 decays. We suppress this jet- 
like background in favor of the spherical signal events, us- 
ing several event-shape variables as in Ref. || . To max- 
imize separation between signal and background, these 
variables are combined in neural networks that are sep- 
arately optimized for each decay mode. Each network 
is trained using Monte Carlo (MC) events, and is val- 
idated on statistically independent MC samples. Cuts 
are made on the neural-network output to suppress con- 
tinuum background. The toes and AA* distributions of 
data are shown in Fig. \I\ior all four A* decay modes. 

The remaining background includes that from BB 
events, which is dominated by B — * A s 7 decays, where 
X s represents hadronic final states other than A* . If one 
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FIG. 1: mES and AE* distributions for the B — > K*^ candidates. The points are data, and the solid and dashed curves show 
the projections of the complete fit and the background component alone, respectively. The fits used to extract the signal yields 
are described in the text. 



TABLE I: The signal efficiency e, the fitted signal yield Ns, the branching fraction B, and the CP-asymmetry A for each decay 
mode. The combined branching fractions and CP-asymmetries for B° — > K*°~f and for B + — > K* + 'y are also shown. Errors 
are statistical and systematic, with the exception of e and Ns, which have only systematic and statistical errors respectively. 
The detailed systematic errors are listed in Table El 



Mode 


e(%) 


N s 


£(xl0" b ) 


Combined ^xlO" 53 ) 


A 


Combined A 




24.4T1.4 


583T30 


3.92±0.20±0.23 


|3.92±0.20±0.24 
|3.87±0.28±0.26 


-0.069 ± 0.046 ±0.011 , 




K a TY° 


15.3T1.9 
17.4T1.6 


62±15 
251T23 


4.02±0.99±0.51 
4.90±0.45±0.46 


0.084±0.075±0.007 


> -0.013 ±0.036 ±0.010 


K a 7V + 


22.1T1.4 


157±16 


3.52±0.35±0.22 


0.061±0.092±0.007 





or more particles escape detection, X s may be incorrectly 
reconstructed as K*, leading to a value of toes near the 
B meson mass, but with AE* distinctly negative. 

For each decay mode, the signal yield and asymme- 
try A (except for the K s tt° mode) are simultaneously ex- 
tracted using an extended unbinned maximum likelihood 
fit, 

3 N 3 

C = exp (- «i) ' [II £ N i V @i 5 » 

i— 1 j— 1 i—1 

to the two-dimensional distribution of toes and AE* 
with three hypotheses (index i): signal, continuum back- 
ground, and B background. The probability density 
function (PDF) Vixj, &i) for each of the three hypothe- 
ses is the product of individual PDFs of the fit variables 
Xj = (mES, AE*). on are the shape parameters for the 
PDFs described below. In the three self-flavor-tagged 
modes (K+tt~, K+tt°, and K s ir+), N. t = \(\-TAi)ni, 



where rii and Ai stand for the total yield and CP- 
asymmetry of signal, continuum background, and B 
background, while in the K s ir° decay mode, Ni = rij. 
The bottom-quark flavor, T , is defined as —1 for b quarks 
and ±1 for b quarks. In the K + ir~ mode, mistagging is 
possible if both the pion and kaon are misidcntified, but 
this probability is negligibly small. We assume that the 
CP asymmetry of the B background and that of the con- 
tinuum background are the same. 

To reduce systematic errors, most of the fit parameters 
for the signal and for the continuum background are de- 
termined by a fit to data. For continuum background, the 
AE* distribution is modeled by a first-order polynomial 
function with the exception of K s ir + , where a second- 
order polynomial is used. The tties distribution for con- 
tinuum background is modeled with an ARGUS function 
[l4| . In the K + tt° decay mode, the continuum back- 
ground shape is simultaneously fit to the off-resonance 
data to obtain a stabler fit. For the B background, the 
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FIG. 2: mj{, spectra for the different decay modes for events 
in the signal region after background subtraction using side- 
bands in tties and AE* . The points are data and solid curves 
represent relativistic p-wave Breit-Wigner line shapes with 
masses and widths of K* taken from Ref . |l3| . 



Gaussian distribution used for AE* and the Novosibirsk 
function |l5| used for ities have all shape parameters 
fixed to values determined from MC. The signal AE* 
distribution is modeled as a Crystal Ball function [lif . 
which is a Gaussian distribution with a low-side power- 
law tail that is fixed using MC. The rriES distribution for 
signal is modeled as a Gaussian function, except for the 
K + tt° decay mode, where a Crystal Ball function, with 
tail parameters fixed using MC fits, is used to accommo- 
date a low-side tail due to the ir° energy lost from the 
calorimeter. The same low-side tail in the K s tt° decay 
mode is ignored due to the small number of events in this 
mode. 

Correlations between m^s & n d AE* distributions 
could introduce a bias in the signal yields. To study this, 
randomly selected events from our detailed MC simula- 
tion of the signal were mixed with background events 
generated using the PDF from the fit. In this way we 
determined that the K + n~ efficiency must be corrected 
by multiplying it by 0.98. For the K^n", K + n° : and 
Kgir + modes, the corresponding numbers are 0.91, 0.96, 
and 0.96. The error in this fit bias due to MC statistics is 
included as a systematic uncertainty. These MC studies 
also indicate that correlations between the B background 
and the continuum background fit yields do not affect the 
fitted signal yield. 

The projections of the maximum likelihood fits on tues 
and AE* are shown in Fig. 2] for each decay mode. Fig- 
ure [21 shows that the background-subtracted Kir invari- 



TABLE II: Fractional systematic uncertainties on the branch- 
ing fractions B and absolute systematic uncertainties on CP 
asymmetry A. 

Systematic errors on B(%) 



Description 


K TV 


7/0 






Number of B events 


1.1 


1.1 


1.1 


1.1 


H 1 


O 1 

ZA 




Z.4 


O 1 

z.4 


1 I'M ( 4 L" 1 II (T oHl^lDH/ 1 !! 

xictcKiiig eiiiciciicy 


1.6 




0.8 


u.o 


V^IldlgeLL JJclI LICie ILLeilLlllCcLLUJIl 


1.0 




1.0 


1 
r .u 


Photon efficiency 


2.5 


7.6 


7.6 


r. 
Z . 


Photon isolation cut 


2.0 


2.0 


2.0 


o n 


/I , // VfclU 


1.0 


1.0 


1.0 


1 

J- .V 


K s efficiency 




3.0 




o.O 


Neural network 


3.0 


3.5 


2.7 


2.8 


PDF parameterization 


2.2 


7.3 


2.7 


1.4 


MC statistics/fit bias 


0.9 


3.2 


2.4 


1.6 


Total 


5.8 


12.3 


9.4 


6.3 




Systematic errors on 


A (%) 


Tracking efficiency 


0.35 




0.25 


0.25 


Charged particle identification 


1.00 




0.55 


0.53 


Nuclear interaction asymmetry 


0.20 




0.35 


0.15 


B-background asymmetry 


0.25 




0.25 


0.25 


Total 


1.1 




0.7 


0.7 



ant mass distributions agree well with the expected K* 
resonance shape. This confirms that the signal is consis- 
tent with coming from only true K* decays. 

Table B shows signal efficiencies, yields from the fits, 
and branching fractions (£>) calculated using our recent 
measurement jl7| of the production ratio of charged and 
neutral B events, R +/0 = r(e+e" -> B + B-)/T(e + e~ -> 
B°B°) = 1.006 ± 0.048 at Vs = M r(4S) . 

Combined values of B{B° K*°j) and B{B+ -> 
K* + j), which are also shown in Table „ are calculated 
taking into account correlated systematic errors between 
modes. We further combined these measurements, us- 
ing the lifetime ratio t b +/t b o = 1.083 ± 0.017 an d 
our measurement of R + /° , to find the isospin asymmetry, 
A _ = 0.050 ± 0.045(stat.) ± 0.028(syst.) ± 0.024(i?+/°), 
which corresponds to an allowed region of —0.046 < 
A _ < 0.146 at the 90% confidence level. We also present 
a combined A measurement in Table„ which corresponds 
to an allowed region of -0.074 < A{B -> K*j) < 0.049 
at the 90% confidence level. 

The systematic error on the branching fraction for each 
mode is shown in Table [n] Most of the uncertainties 
are determined as in our previous analysis 6], so we 
provide details only for the new procedures used. The 
neural-network inputs are generally independent of the 
fully reconstructed B — * K*j candidate, so we determine 
their efficiencies and systematic uncertainties with high- 
purity control samples with reconstructed B~ — > D°ir~ 
and B° — > D~ir + . The "PDF parameterization" error 
comes from MC studies of our fitting procedure, in which 
we estimate the uncertainty incurred by fixing parame- 
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ters in the continuum and B background models. This 
includes uncertainty in the inclusive branching fraction 
and spectral shape of B — > Xgj. 

The systematic uncertainties in the measurement of 
A are also shown in Table |n] The first three contri- 
butions arise from potential particle-antiparticle asym- 
metries in the detector response, including differences 
in interaction cross-sections for K + and K~, and for 
7r + and tt~ (estimated with a method similar to that 
used in Ref. ^3)- The uncertainty due to a possible 
asymmetry in the B background, which is dominated by 
B — > X s j, is estimated by varying our recent measure- 
ment of A(B -> X sl ) 0* within its errors. 

We conclude that both the isospin- and CP- 
asymmetries in B — > K*-f decay processes are consis- 
tent with SM predictions. The branching fractions mea- 
sured are also consistent with SM-based calculations and 
are more precise than those predictions. These measure- 
ments are consistent with previous results 0, 0> 0] • 
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